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ABSTRACT 

We report high sensitivity sub-arcsecond angular resolution observations of 
the massive star forming region DR21(0H) at 3.6, 1.3, and 0.7cm obtained with 
the Very Large Array. In addition, we conducted observations of CH3OH 44 GHz 
masers. We detected more than 30 new maser components in the DR21(0H) 
region. Most of the masers appear to trace a sequence of bow-shocks in a bipolar 
outflow. The cm continuum observations reveal a cluster of radio sources; the 
strongest emission is found toward the molecular core MMl. The radio sources in 
MMl are located about 5" north of the symmetry center of the CH3OH outflow, 
and therefore, they are unlikely to be associated with the outflow. Instead, the 
driving source of the outflow is likely located in the MM2 core. Although based 
on circumstantial evidence, the radio continuum from MMl appears to trace free- 
free emission from shock-ionized gas in a jet. The orientation of the putative jet in 
MMl is approximately parallel to the CH3OH outflow and almost perpendicular 
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to the large scale molecular filament that connects DR21 and DR21(0H). This 
suggests that the (accretion) disks associated with the outflows/jets in the DR21 

- DR21(0H) region have symmetry axes mostly perpendicular to the filament. 

Subject headings: HII regions — ISM: molecules — masers — stars: formation 

— ISM: individual (DR21(0H)) 



1. Introduction 

Located about 3' north of DR21, the DR21(0H) region, also known as W75S, is a site 
of dense molecular gas within the Cygnus-X complex. Both DR21 and DR21(0H) have been 
extensively studied in the infrared, sub-millimeter, and millimeter bands. Both are regions 
of massive star formation, but DR21 is in a more evolved state (as evident from the presence 
of ultracompact H II regions) while DR21(0H) is in an earher phase, in which the massive 
(proto)stars have not yet substantially ionized the surrounding molecular gas. 

The DR21(0H) region has four principal condensations, i.e., DR21(0H)-N, M, W, and 
S (Mangum et al. 1992). DR21(0H)-M contains two main clumps (MMl and MM2) that 
are warm, 60 and 30 K, and massive, 350 and 570 Mq, respectively (Mangum et al. 1992; 
see also Liechti & Walmsley 1997). The total mass and luminosity of the region are about 
3 X 10^ Mq and 5 x 10"^ Lq (Chandler et al. 1993a, 1993b). Star formation activity is indicated 
by the presence of millimeter continuum sources (e.g., Mangum et al. 1991; Chandler et al. 
1993a), a high velocity molecular outflow detected in CS J=5-4 (Richardson et al. 1994) 
and various species of masers (e.g., Plambeck & Menten 1990; Kurtz et al. 2004; Fish et 
al. 2005). Particularly interesting is the outflow traced by CH3OH 44 GHz masers: sub- 
arcsecond angular resolution observations by Kogan & Slysh (1998) show the presence of 
several CH3OH 44 GHz masers grouped in a bipolar structure. The morphology of the 
blueshifted masers resembles a bow shock. Higher sensitivity CH3OH 44 GHz observations 
by Kurtz et al. (2004) confirm the Kogan & Slysh (1998) results and reveal a possible second 
bow shock. 

In summary, numerous observations of DR21(0H) clearly identify it as a very young 
massive star forming region, with multiple OB-star clusters forming. All of the clusters 
appear to be in a pre-ultracompact H II region stage. The unambiguous evidence of massive 
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star formation and the relatively nearby distance (~ 2kpc; Odenwald & Schwartz 1993) 
make DR21(0H) an ideal laboratory to study early stages of massive star formation. 

Extensive work by several groups (e.g., see Eisloffel et al. 2000) has shown that star 
formation regions usually contain numerous weak centimeter continuum sources. These 
sources trace a variety of astrophysical phenomena, including thermal jets, stellar winds, 
and gyro-synchrotron emission from young stellar objects. The location and characterization 
of these weak radio sources has proved extremely valuable to our understanding of star 
formation phenomena (e.g., Zapata et al. 2004; Ignace & Churchwell 2004; Girart et al. 
2002; Reipurth et al. 2002). 

Despite a wealth of molecular line and millimeter continuum observations of DR21(0H), 
it has never been reliably imaged with high sensitivity and angular resolution in the centime- 
ter continuum. Relatively low quality VLaI^ archive continuum observations were reported 
by Argon et al. (2000) and Fish et al. (2005), and the most sensitive cm continuum data 
available in the literature has 3a sensitivities of 10 - 15 mJybeam"^ (Johnston et al. 1984; 
Mangum et al. 1992). In this work we present high angular resolution {Ogyn < 1") and high 
sensitivity {rms < ImJy) observations of the radio continuum at 3.6, 1.3, and 0.7 cm, as 
well as observations of the CH3OH 44 GHz masers in DR21(0H). 



2. Observations and Data Reduction 

We used the VLA to conduct radio continuum observations of DR21(0H) at X (3.6 cm), 
K (1.3 cm) and Q (7 mm) bands. In addition, we report observations of CH3OH 44 GHz 
masers. Details of the observations are summarized in Table l|f|. 

The 3.6 cm and 1.3 cm observations were conducted using the standard VLA continuum 
mode, i.e., 50 MHz bandwidth, 4IF mode. The quasar 3C48 was used as primary (flux) cali- 
brator and J2007+404 was used as secondary (phase) calibrator. The calibration cycle of the 
3.6 cm and 1.3 cm observations was ~ 600/60 and 220/40 s (source/calibrator), respectively. 
All data reduction was done in AIPS following standard procedures. 

The 7 mm continuum was observed simultaneously with the CH3OH masers using two 
IF pairs. The goal of this setup was to cross-calibrate the radio continuum observations 



^The Very Large Array (VLA) is operated by the National Radio Astronomy Observatory (NRAO), a 
facility of the National Science Foundation operated under cooperative agreement by Associated Universities, 
Inc. 

•^The observations reported in this paper are from the NRAO proposals AH745 and AK586. 
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using the CH3OH masers to minimize phase decorrelation due to tropospheric variations. 
One pair of IFs was set to a narrow bandwidth of 3.125 MHz (63 channels, 48.83 kHz) to 
observe the CH3OH uq = 44.06943 GHz (7o ~ 61, A"*") hne; the second pair of IFs was used 
with a bandwidth of 25.0 MHz to observe the radio continuum (excluding the maser lines). 
The narrow band IF was centered on Vlsr = 0.0 km s~^. The CH3OH spectra were Hanning 
smoothed to a final channel width of 97.6 kHz (0.66 km s~^). We observed 3C286 as primary 
(flux) calibrator and 1229+020 as bandpass calibration; in addition, J2007+404 was observed 
approximately every hour to correct for phase offsets between the CH3OH and continuum 
passbands. Assuming the position reported in Kurtz et al. (2004) for the brightest CH3OH 
maser in the region, we self-calibrated the peak maser channel and transferred the calibration 
table to all CH3OH channels, as well as to the broadband (25 MHz) continuum IFs. All 
calibration and imaging was done using the NRAO package AIPS. 



3. Results &: Discussion 

We detected radio continuum emission in all three bands (3.6 cm, 1.3 cm, and 7 mm) as 
well as CH3OH 44 GHz masers. Figure 1 shows the radio continuum contour maps; Table 2 
lists the parameters of the final continuum images and in Table 3 we report the observed 
properties of all radio continuum sources. Figure 2 shows the combined CH3OH spectrum 
(measured from a 35"x 11" box centered between the MMl continuum sources and 3.6cm 
source Rl), and Figure 3 shows the integrated intensity (zero velocity moment) image of 
the CH3OH emission. In Figure 4 we show the velocity field (first velocity moment) of the 
CH3OH maser arcs seen in Figure 3. To facilitate comparison between the figures, we also 
show in Figure 4 the 3.6 cm radio continuum emission (contours; see Figure 1). CH3OH line 
parameters are listed in Table 4. 

3.1. CH3OH 44 GHz Masers: Outflow in DR21 (OH) 

The CH3OH masers were observed in this work with a two-fold goal: to cross-calibrate 
the 7 mm radio continuum and to observe the masers with higher angular resolution {Ogyn ~ 
0.6" versus 1.7") and sensitivity (rms ~ 5mJybeam~^ versus 50mJybeam~^) than the ob- 
servations reported by Kurtz et al. (2004). Our observations however, have poorer spectral 
resolution than those of Kurtz et al. (2004), i.e., 0.7 versus 0.17km s~^. Given the difference 
in spectral and angular resolution, as well as flux density calibration uncertainties (see Kurtz 
et al. 2004), we consider that our data are not suitable for a variabihty study of the masers 
by comparing these data sets. 
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Figure 3 shows the integrated intensity distribution of CH3OH masers in DR21(0H). 
We detected more than 30 new maser spots. Figure 4 shows the velocity distribution. As 
noted by Kurtz et al. (2004) and Kogan & Slysh (1998), the CH3OH 44 GHz masers trace 
a bipolar structure, with blueshifted masers located predominately in the eastern lobe, and 
red shifted masers located in the western lobe. As proposed by Kurtz et al. (2004, see 
also Plambeck & Menten 1990 for the case of 95 GHz CH3OH masers), this organized maser 
distribution (in space and velocity) is likely tracing an outflow in DR21(0H). Given the 
mass and temperature of the MM2 core (Mangum et al. 1991), the presence of an east- 
west molecular outflow (e.g., Lai et al. 2003; Chandler et al. 1993b) and the widespread 
association of CH3OH Class I masers with young massive stars (e.g., Ellingsen 2006; Pratap 
et al. 2008), the driving source is likely a young massive stellar object somewhere in MM2, 
located toward the center of the bipolar maser distribution. 

The distribution of masers within each lobe is quite intriguing. Kogan & Slysh (1998) 
found that the CH3OH 44 GHz masers in the blueshifted group are distributed along an arc, 
suggestive of a bow-shock origin. The higher sensitivity observations of Kurtz et al. (2004) 
confirm the morphology in the blue lobe and show the presence of a possible second arc. Our 
still higher sensitivity observations enable a more detailed characterization of the CH3OH 
maser distribution in DR21(0H). 

Our data confirm the two-arc maser distribution of the eastern (blue) lobe and reveal 
that the masers in the western lobe are also grouped in two arc-like structures. As explicitly 
suggested by thin lines in Figure 4, the CH3OH masers appear to trace two nested bow-shocks 
in the outflow. This distribution is remarkably similar to the jet bow-shock outflow model 
shown by Arce et al. (2007). The occurrence of sequential bow-shocks (outflow events) is 
also reminiscent of the morphology of the outflow in DR21 (Smith et al. 2006). Clearly, the 
sky distribution of the CH3OH masers traces a 2D projection of a 3D structure; the masers 
might originate in a cut through parabolic-like outflow lobes where the velocity coherence 
and beaming are favorable for maser amplification toward our line of sig htQ 

Given the symmetric bow-shock morphology and the low velocity gradient between the 
red and blueshifted maser groups (< 10km s~^. Figure 4), it appears that the outflow is 
almost in the plane of the sky. Assuming that the CH3OH masers are tracing two individual 
(nested) outflow events, and assuming, for simplicity, a constant expansion velocity of 100 km 
s~^ for each lobe (e.g., as observed in H2O masers in IRAS 20126-1-4104, Moscadelli et al. 
2005), then the dynamical ages are ~ 900 yr and ~ 1500 yr for the inner and outer shocks, 
respectively (a distance of 2kpc was assumed), i.e., the lower limit of the outflow age is 
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~ 1500 yr. Thus, the outflow could be quite young — younger than the 'pulsed' outflow 
in Cep-A (e.g., Bally 2008). However, we note that high angular resolution (VLBI) proper 
motion studies are required to determine the velocity of the CH3OH masers in the region. In 
particular, CH3OH molecules may not survive in high velocity J-shock environments (e.g., 
Garay et al. 2002). Alternatively, the CH3OH masers may be tracing a slower shock traveling 
along the high speed flow and/or the working surface between two shocks]^ 

A complication to the single-outflow interpretation of the CH3OH masers in DR21(0H) 
is the misalignment between the arcs, particularly with respect to the eastern arc (Figure 4; 
i.e., a change in position angle of ~15° east of north). Such a change in position angle 
could result from jet precession combined with inhomogeneities in the molecular gas that 
is interacting with the flow, similar to the large scale outflow in G192. 16—3.82 (Devine et 
al. 1999). Nevertheless, despite the misalignment, the outflow seems to be quite collimated, 
with an apparently constant opening cross-section at scales greater than 10^ AU (~0.05pc) 
from the symmetry center. 



3.2. Radio Continuum: A Cluster of Radio Sources 

In order to investigate the driving source of the CH3OH maser outflow in DR21(0H), 
we conducted sub-arcsecond radio continuum observations at 3.6 cm, 1.3 cm, and 7mm (Fig- 
ure 1). We detected eight radio continuum sources (see Table 3). In Figure 4 we show the 
CH3OH velocity field superimposed on the 3.6 cm continuum map. One of the 3.6 cm con- 
tinuum sources (Rl) is located close to the projected center of the CH3OH maser outflow; 
it may mark the position of the massive (proto)star that drives the outflow. 

Alternatively, DR21(0II)-R1 could mark the position of a young low-mass star. The 
radio emission could be gyro-synchrotron radiation from an active magnetosphere, just as 
the low mass companion of 9^ Ori A (Felli et al. 1993) or the variable radio source in 
IRAS 20126-1-4104 (Hofner et al. 2007). This hypothesis is supported by the detection of 
other weak radio sources (R5, R6 and possibly R2; Figure 1), that together with Rl may 
represent the high end of the (radio) luminosity function of a cluster of low mass objects 
forming among the young massive stars. 

In Table 5 we list the spectral index between 3.6 and 1.3 cm of all sources detected at 
both wavelengths. The spectral indexes of Rl, R2, R5, and R6 (a different interpretation 
of the radio emission from R3 and R4 is given below) are between —0.6 and 0.4, which 
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is similar to the range of spectral index variability of the low mass star in 9^ Ori A (FeUi 
et al. 1993). Further observations of the weak population of radio sources in DR21(0H) 
are necessary to assess whether the radio emission is variable as expected from the gyro- 
synchrotron hypothesis, and with this, to study how low-mass pre-main-sequence stars (class 
I/ll) and very young class high-mass objects can coexist considering their clearly different 
evolutionary timcscales. 

In contrast to the radio emission from Rl, R5, and R6 (and possibly R2), the radio 
sources MMl-NW, MMl-SE, R3 and R4 appear to be related to the massive young stellar 
object in MMl. Wc propose that these radio sources arc tracing shock- ionized gas in the jet 
from a young massive star. Wc prefer this interpretation because: 

1. MMl-NW, SE, R3 and R4 are coUinear (southeast-northwest orientation), which suggests 
a common driving source. 

2. In principle the similar flux densities of MMl-NW and SE could be explained by radio 
emission from two hypercompact H II regions that are in a very similar evolutionary state. 
However, the morphology of the two radio sources is also similar; after deconvolution, both 
sources appear to be elongated at approximately the same position angle (~ 40°; see Table 
3). This position angle is almost perpendicular to the MMl-NW to R4 orientation (i.e., 
orthogonal to the putative jet; see Figure 1), and may arise from a bow shock aligned with 
the MMl-NW to R4 axis. The elongation would result from our edge-on view of the shocks, 
corresponding to the broad tail of a relatively large-opening-angle bow shock, emitting free- 
free radiation. Given the symmetry, we consider it more hkely that MMl-NW and SE are 
manifestations of a single radio jet. 

3. The spectral indexes between X and K bands are very similar for MMl-NW and SE 
{a ~ 0.8); in addition, the spectral index of R4 is also similar to those of MMl-NW and SE 
within the errors (Table 5). A spectral index of ~ 0.8 is consistent with thermal emission 
from finite stratified density gas being ionized by UV photons from the cooling region of a 
shock (e.g., Ghavamian & Hartigan 1998). 

The weak radio emission supports the interpretation that the most massive stellar object 
in MMl is a high mass protostellar object (HMPO; e.g., Beuther et al. 2007). 

In Figure 1 we show the location of 6.7 GHz CH3OH masers (Harvey-Smith et al. 2008) 
and H2O masers (Palmer & Goss, in prep.). The H2O masers are distributed throughout the 
region, and thus are unlikely to be driven by a single young stellar object. Two of the weak 
radio sources (Rl and R6) may be associated with H2O masers, and another H2O maser is 
found almost at the center of MMl-NW/SE. The detection of a H2O maser at the center of 
the double radio source is reminiscent of IRAS 20126-1-4104 (Hofner et al. 2007; Moscadelli 
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et al. 2005) where a clump of H2O masers is also located toward the center of a double radio 
source. The radio continuum emission associated with the H2O masers in IRAS 20126+4104 
is likely due to shock-ionized gas in the jet from a young massive star (see also §3.5). 

3.3. The Nature of the 7 mm Emission 

As evident from Figure 1, the 7 mm source has the elongation expected from superpo- 
sition of emission from MMl-NW, SE, R2, R3, and R4 when observed with lower angular 
resolution. The southeast - northwest elongation is also consistent with the 2 cm continuum 
map shown by Fish et al. (2005). To further investigate the nature of the 7 mm source 
we convolved the 3.6 and 1.3 cm observations to the synthesized beam of the 7 mm map. 
In Figure 5 we show the (radio) flux density distribution of MMl after convolution, i.e., 
MMl-NW + SE + R2 + R3 + R4. The 7 mm flux density point lies slightly above the 
power law fit. This indicates that at 7mm the flux is dominated by free- free emission with 
a small (~ 1.4 mJy or 10% of the total flux density) excess, most likely due to thermal dust 
emission. 



3.4. Infrared Counterparts 

The DR21/W75 region (including DR21(0II)) has been extensively studied in the in- 
frared (e.g., Davis et al. 2007; Kumar et al. 2007; Smith et al. 2006; Marston et al. 2004). 
In Figure 6 we show the Spitzer IRAC image of the region (3.6 /im blue; 4.5 fim green; 8.0 //m 
red)lj. These IRAC observations have been discussed in a number of papers (e.g., Kumar 
et al. 2007) and are presented here to facilitate comparison between the CII3OH, radio 
continuum and infrared data. 

Figure 6 shows a number of sources with 4.5 /im excess particularly along the large scale 
north-south (dark) filament (see Kumar et al. [2007] for a larger scale view of the region). 
Excess in the 4.5 /um band is a known tracer of shocks (e.g.. Smith et al. 2006; Cyganowski et 
al. 2008), thus. Figure 6 exemplifies that star formation in DR21(0II) is not limited to the 
mm cores but is found throughout the region. The eastern CH3OH arcs may be associated 
with 4.5 fim excess emission, however there is no clear evidence for 4.5 /xm excess toward the 
western arcs. 



^Data were retrieved from the 5piteer archive (http://irsa.ipac.caltech.edu/applications/Spitzer/Spitzer/ ) 
project PID-1021, "DR21 and its Molecular Outflow", Marston et al. (2004). 
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Shocked gas in the region has also been studied by Davis et al. (2007) based on H2 
2.12 fim observations. These authors did not find H2 emission corresponding to the CH3OH 
masers. This might be due to the high optical depth at 2.12 /xm for these deeply embedded 
regions in the vicinity of MMl and MM2. However, slightly outside the area shown in 
Figure 4, Davis et al. (2007) found two H2 regions that are somewhat aligned with the 
maser outfiow (B7-1 and B5-1; see Table 3 and Fig. A. 2 of Davis et al. [2007]). 

The MMl radio continuum sources are almost coincident with the peak of an extended 
"red" object (see Figure 6); the infrared counterpart of MMl was also detected at 24 ^m 
(saturated) and 70 fim with MIPS. Figure 6 also shows a "blue" star located very close to 
the symmetry center of the CH3OH outfiow. The star is clearly detected at 3.6 and 4.5 ^m 
(brighter at 3.6 fim) and tentatively detected at 5.8 //m blended with the MMl emission. 
The star is also detected in J, H, and K bands (2MASS), and in optical DSS images. 

The near-infrared (NIR) colors of the "blue" object are intriguing. The 2MASS colors 
are J — H = 0.32 and H — Kg = 0.11, which, if placed into a NIR color-color diagram (e.g., 
Hillenbrand et al. 1998, Meyer et al. 1997, Marston et al. 2004, Apai et al. 2005) are 
consistent with an early F-type main sequence star with ~1 mag of optical extinctioiu. If 
it is a FO star, then to attain V = 13 mag (as roughly observed), such a star would be at a 
distance of ~700 pc, and hence would not be associated with the DR21 region. However, the 
Spitzer colors are inconsistent with this interpretation. If the optical extinction were 1 mag, 
a Spitzer IRAC color of [3.6]-[4.5] = 0.53 could not be caused by extinction, because this 
color indicates much higher extinction values [Ay >10 mag). Hence, the blue object could 
still be associated with the DR21 region, and other circumstances (unresolved multiplicity, 
existence of a transitional circumstellar disk, etc.) might govern the infrared colors. Based 
on the current data, it is hard to decide whether the blue object is the driving source for the 
CH3OH outfiow. An in-depth investigation of the nature of this near-IR source is beyond 
the scope of this work, and should be the topic of a future study. 



3.5. A New Observational Phase during the Formation of Massive Stars 

It is interesting to compare the characteristics of the radio emission in DR21(0H) with 
that of other young massive stellar objects, in particular G31.41+0.31 and IRAS 20126+4104 
(Araya et al. 2003, 2008; Hofner et al. 2007). In all three cases, there are multiple radio 



^More specifically, one magnitude of optical extinction corresponds to extinctions of 0.28 mag in J band, 
0.18 mag in H band, and 0.11 mag in K band (Mathis 1990). Correcting the 2MASS colors for these 
extinction values and based on Ducati et al (2001), the colors are consistent with those of an FO - F2 star. 
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continuum sources toward the center of molecular cores. Some of the continuum sources are 
consistent with radio emission from shock-ionized gas and not due to direct photo-ionization 
from an embedded massive star (see for example figure 1 of Hofner et al. 2007). The 
intermediate mass object GGD27 (HH80/81) is a clear example of ionized knots moving 
along a collimated jet as shown by proper motion studies (Marti et al. 1995, 1998). 

We propose that there is a period during the formation of massive stars when (possibly 
recurrent) jet events interact with the molecular core and ionize the medium via shocks. 
Such a process could produce the double radio sources we have detected in the core of 
several young massive star forming regions. The ionization process would eventually become 
dominated by the development of hypercompact and then ultracompact H II regions. 

3.6. Further Implications 

The data reported in this paper show regularities in the orientation of outflows among 
themselves and also with respect to the large scale distribution of dust and gas in the region. 
The angular distance between the MMl radio emission and the geometrical center of the 
CH3OH outflow (see Figure 4) is too large to assume that the driving source of the outflow 
is located in MMl. Thus, our data indicate the presence of two distinct young massive stellar 
objects (or systems) with nearly parallel outflows/jets; i.e., both the putative MMl-NW to 
R4 jet and the CH3OH outflow have similar position angles (~ 130°). This suggests that 
the orientation of jets/outflows is. to flrst order, determined by the large scale properties of 
the parent cloud. This conclusion is further supported by the powerful outflow from DR21 
which also has a dominant east- west orientation (e.g., Marston et al. 2004). The statistically 
signiflcant orientation of the outflows in an east-west direction was also discussed by Davis 
et al. (2007). Other examples of quasi-parallel outflow from different young massive stellar 
objects within the same star forming region include Cep-A and IRAS 05358-1-3543 (Beuther 
et al. 2002; Bally 2008). 

The complete DR21 to DR21(0H) region is part of a long molecular cloud fllament (or 
sheet seen edge-on) oriented north-south (e.g., Marston et al. 2004; Motte et al. 2005, 2007). 
The outflows are mostly perpendicular to the large scale north-south filament, and approx- 
imately aligned with the large scale direction of the magnetic field through the filament, 
which, as discussed by Vallee & Fiege (2006), shows a predominant east- west orientation 
(see also Lai et al. 2003). If the outflows arise from (proto)stars with accretion disks, then 
the disks would tend to be oriented edge-on along the filament, with disk symmetry axes 



oriented east-west 



4. Summary 

We report sub-arcsecond observations of the radio continuum at 3.6, 1.3, and 0.7 cm, 
and also of the CH3OH 44 GHz masers in the DR21(0H) region. We detected a cluster of 
radio sources; the two brightest sources are located toward the molecular core MMl and 
appear to trace radio emission from shock-ionized gas. 

We detected more than 30 new CH3OH 44 GHz maser spots in DR21(0H). Our obser- 
vations confirm and delineate a bow-shock distribution of the masers. The masers show a 
well-separated (in space and velocity) bipolar outflow almost symmetrically located with re- 
spect to the MM2 core, which is likely to harbor the driving source of the flow. The CH3OH 
masers appear to trace two different outburst events. The sequential bow-shock morphology 
observed in DR21(0H) is reminiscent of that of the DR21 outflow (Smith et al. 2006). 

The driving source of the CH3OH outflow in DR21(0H) is unlikely to be the one re- 
sponsible for the radio emission in MMl. More likely, two different young massive stellar 
objects are driving outflows/jets at very similar position angles. At even greater scales, the 
prominent outflow from DR21 also has an ~east-west orientation. As also noted by Davis et 
al. (2007), throughout the large scale filament in the DR21-DR21(OH) region we find that 
the outflows/jets are quite parallel in a ~ east-west direction, which is perpendicular to the 
(north-south) filament. This suggests that the accretion disks form mostly edge-on along 
the large scale filament, i.e., disk symmetry axis orthogonal to the filament elongation. 

E.A. acknowledges support from an University of New Mexico postdoctoral fellowship 
and a National Radio Astronomy Observatory (NRAO) Jansky fellowship. We thank L. 
Deharveng for a discussion regarding the IR emission in DR21(0H) and an anonymous referee 
for a thorough review and comments that improved the manuscript. We also thank P. Palmer 
and M. Goss for providing the locations of the H2O masers. This research made use of the 
NASA's Astrophysics Data System, and archival data from the Spitzer Space Telescope and 
the Digitized Sky Survey. The Digitized Sky Surveys were produced at the Space Telescope 
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potential counterexample is EROS in the DR21(0H)N region, where Harvey-Smith et al. (2008) and 
Harvey-Smith & Soria-Ruiz (2008) report a possible Keplerian disk traced by CH3OH masers and oriented 
east-west. However, we consider that future studies may be needed to confirm their result. 
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Table 1. VLA Observations 



Parameter 


X-Band (3.6 cm) 


K-Band (1.3 cm) 


Q-Band (7 mm)" 


CH3OH 44 GHz" 


Date 


Nov 16, 2004 


May 05, 2005 


Sep 07, 2001 


Sep 07, 2001 


VLA Configuration 


VLA-A 


VLA-B 


VLA-C 


VLA-C 




20 39 01.000 


20 39 01.000 


20 39 00.900 


20 39 00.900 


Dec'' 


42 22 40.00 


42 22 40.00 


42 22 47.00 


42 22 47.00 


uo (GHz) 


8.46 


22.4 


44.1 


44.06943 


Correlator Mode 


4IF Cont. 


4IF Cont. 


4IF Line 


4IF Line 


Bandwidth per IF (MHz) 


50 MHz 


50 MHz 


25.0 (2 Pol) 


3.125 (2 Pol) 


Flux Density Calib. 


3C48 


3C48 


3C286 


3C286 


Assumed S,/ (Jy) 


3.15 


1.12 


1.43 


1.43 


Phase Calib. 


J2007+404 


J2007+404 


J2007+404= 


J2007-|-404= 


Measured Si/ (Jy) 


2.54±0.01 


1.78d=0.02 


1.34±0.04 


1.35±0.04 



"The Q-Band (7 cm) and CH3OH 44 GHz observations were conducted simultaneously using a 4IF mode. 
Self-calibration solutions from the CH3OH maser observations were used to calibrate the Q-Band continuum. 

''Phase tracking center (J2000). 



"^Used to correct phase offsets between the wideband and CH3OH spectral line band. 
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Table 2. Parameters of the Radio Continuum Images 



Parameter 


X-Band {3.6 cm) 


K-Band (1.3 cm) 


Q-Band (7 mm) 


Syn. Beam 


0'.'20 X Of.' 19 


0'.'25 X Of.' 25 


Of.' 92 X Of.' 74 


Syn. Beam P.A. 


-18.2° 


62.1° 


-74.9° 


rms (;iJyb~^) 


19 


36 


400 
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Table 3. Radio Continuum Sources" 



Source 


Band 


RA(J2000) 


Dec(J2000) 


Size^ RA. 










20''39™ 


42°22' 


(AU,°) 


(mjy beam~l) 


(mJy) 


DR21(OH)-Rl 


3.6 cm 


00':612{1) 


43'.' 44(1) 


440x100, 86° 


0.17(2) 


0.26(4) 




1.3 cm 


00! 606(4) 


43('45(3) 


670x <160, 59° 


0.14(4) 


0.20(8) 


MMl-NW 


3.6 cm 


00=9763(1) 


48"980(2) 


290 X <50, 34° 


1.14(2) 


1.40(4) 




1.3 cm 


00?9765(1) 


48'.'981(1) 


200x80, 37° 


2.70(4) 


2.95(7) 


MMl-SE 


3.6 cm 


01?0079(1) 


48'.'707(2) 


280x155, 46° 


1.18(2) 


1.55(4) 




1.3 cm 


01!0079(1) 


48'.'695(1) 


235x225, 51° 


2.83(4) 


3.43(7) 


MM1= 


7 mm 


01?02'* 


48'.'8'^ 


1830 X <460, 123° 


9.4(4) 


13.4(8) 


DR21(OH)-R2 


3.6 cm 


01?042(1) 


48'.'93(1) 


330 X <100, 65° 


0.16(2) 


0.18(4) 




1.3 cm 


01?045(3) 


48'.' 9 1(2) 


440 X <300, 70° 


0.16(4) 


0.20(7) 


DR21{OH)-R3 


1.3 cm 


01?054(2) 


48('43(2) 


690x <300, 111° 


0.19(3) 


0.32(9) 


DR21(OH)-R4 


3.6 cm 


01!083(2) 


48'.'35(2) 


465x420, 94° 


0.14(2) 


0.31(6) 




1.3 cm 


01!077(3) 


48'.'32(2) 


745x315, 97° 


0.21(3) 


0.45(10) 


DR21(OH)-R5 


3.6 cm 


01?303(1) 


49'.'019(8) 


285 X <135, 77° 


0.18(2) 


0.21(4) 




1.3 cm 


01!308(2) 


49'.'05(2) 


< 230x < 90, 144° 


0.16(4) 


0.12(5) 


DR21(OH)-R6 


3.6 cm 


01':337(1) 


52('00(2) 


190x <310, 60° 


0.11(2) 


0.12(4) 




1.3 cm 


01!343(2) 


51'.'98(2) 


180x < 100, 125° 


0.20(4) 


0.18(6) 



''Parameters obtained from a 2D Gaussian fit of the brightness distribution using the task JMFIT in AIPS. 
The format of the error notation is, e.g., xx.xxx(y) = xx.xxx ± O.OOy. 

''Nominal deconvolved size (major and minor axes, and position angle). A distance of 2kpc was assumed 

(Odenwald & Schwartz 1993). We list a quarter of the synthesized beam (0''05~100 AU at X-Band) as upper 
limit of the size in the case of unresolved sources for which JMFIT does not report a limit. 

"^The sources MMl-NW, MMl-SE, DR21(OH)-R2, R3, and R4 are blended at 7 mm. 

''For the self-calibration, we assumed the position of the brightest CH3OH maser as given by Kurtz et al. 
(2004). Thus, the position error of our Q-band observations is dominated by the Kurtz et al. (2004) astrometric 
accuracy of ~0'.'6. 
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Table 4. CH3OH Line Parameters 



Maser ID 


RA (J2000) 


Dec (J2000) 


IS 

(Jyb-1) 


(km s-i) 


AV 
(km 


'a 
3a 

s-i) 


/ Ii,dv 
(Jyb-i kms-i) 


Notes'' 


CH3OH-I 


20 


38 


58 


.68 


42 


22 


22, 


.1 


0.40 


-1.7 


2, 


.0 


0.41 


N 


CH3OH-2 


20 


38 


58 


.92 


42 


22 


27, 


.2 


0.39 


-6.3 


3, 


.3 


0.92 


N 


CH3OH-3 


20 


38 


59 


.29 


42 


22 


48, 


.8 


76.8 


0.3 


3, 


.3 


79.3 


Y-1 


CH3OH-4 


20 


38 


59, 


.32 


42 


22 


49, 


.6 


2.95 


1.0 


2, 


.0 


3.28 


N 


CH3OH-5 


20 


38 


59, 


.35 


42 


22 


47. 


.4 


2.05 


-1.7 


2, 


.6 


2.68 


Y-2 


CH3OH-6 


20 


38 


59, 


.56 


42 


23 


05. 


.3 


0.49 


-1.0 


2, 


.0 


0.57 


Y-3 


CH3OH-7 


20 


38 


59, 


.64 


42 


22 


47, 


.2 


0.092 


-0.3 


2, 


.6 


0.16 


N 


CH3OH-8 


20 


38 


59 


.72 


42 


22 


49, 


.6 


0.27 


-0.3 


2, 


.6 


0.29 


N 


CH3OH-9 


20 


38 


59 


.73 


42 


23 


16, 


.2 


0.27 


-0.3 


2, 


.6 


0.49 


N 


CH3OH-IO 


20 


38 


59 


.73 


42 


22 


47, 


.6 


0.12 


-0.3 


1 


.3 


0.12 


N 


CH3OH-II 


20 


38 


59, 


.76 


42 


22 


45, 


.5 


0.23 


-0.3 


2, 


.6 


0.37 


N 


CH3OH-I2 


20 


38 


59, 


.83 


42 


22 


50. 


.1 


0.14 


-1.7 


3, 


.3 


0.23 


N, 


CH3OH-I3 


20 


38 


59, 


.89 


42 


22 


46. 


.0 


3.01 


-0.3 


3, 


.3 


3.93 


Y-4, 


CH3OH-I4 


20 


38 


59, 


.92 


42 


22 


47, 


.9 


0.12 


-0.3 


3 


.3 


0.18 


N 


CH3OH-I5 


20 


38 


59 


.92 


42 


22 


50, 


.3 


0.21 


-1.0 


3, 


.3 


0.28 


N, 


CH3OH-I6 


20 


38 


59, 


.93 


42 


22 


45, 


.0 


5.27 


0.3 


2, 


.6 


6.07 


Y-6 


CH3OH-I7 


20 


38 


59, 


.99 


42 


22 


35. 


.2 


0.16 


-2.3 


1, 


.3 


0.14 


N 


CH3OH-I8 


20 


39 


00 


.08 


42 


22 


47. 


.1 


0.26 


-1.7 


2, 


.0 


0.38 


N 


CH3OH-I9 


20 


39 


00 


.09 


42 


22 


44. 


.7 


0.14 


-8.3 


1 


.3 


0.16 


N 


CH3OH-2O 


20 


39 


00 


.12 


42 


22 


43. 


.7 


0.38 


-5.6 


2, 


.0 


0.42 


N 


CH3OH-2I 


20 


39 


00 


.13 


42 


22 


47, 


.1 


0.43 


-0.3 


2, 


.6 


0.83 


N 


CH3OH-22 


20 


39 


00 


.14 


42 


22 


48, 


.8 


0.081 


-3.0 


3, 


.3 


0.19 


N 


CH3OH-23 


20 


39 


00, 


.19 


42 


22 


47. 


.7 


3.26 


-1.7 


2, 


.0 


3.48 


Y-7 


CH3OH-24 


20 


39 


00, 


.26 


42 


22 


46. 


.0 


21.3 


-0.3 


3, 


.3 


25.7 


Y-8 


CH3OH-25 


20 


39 


00, 


.29 


42 


22 


47. 


.4 


5.74 


-0.3 


3, 


.9 


7.10 


N 


CH3OH-26 


20 


39 


00 


.33 


42 


22 


48, 


.5 


0.15 


-1.0 


2, 


.0 


0.23 


N, 


CH3OH-27 


20 


39 


00 


.49 


42 


22 


46, 


.3 


0.097 


-3.0 


2, 


.0 


0.14 


N 


CH3OH-28 


20 


39 


00, 


.55 


42 


22 


47, 


.7 


1.08 


1.7 


2, 


.6 


1.47 


Y-9 


CH3OH-29 


20 


39 


00, 


.60 


42 


22 


45. 


.0 


0.12 


-2.3 


2, 


.0 


0.20 


N 


CH3OH-3O 


20 


39 


01, 


.05 


42 


22 


41. 


.6 


0.29 


—1.7 


3, 


.3 


0.40 


N 


CH3OH-3I 


20 


39 


01, 


.05 


42 


22 


17. 


.6 


0.34 


-3.7 


2, 


.0 


0.37 


N 


CH3OH-32 


20 


39 


01, 


.09 


42 


22 


01. 


.0 


2.21 


-2.3 


3. 


.3 


2.51 


N 


CH3OH-33 


20 


39 


01, 


.17 


42 


22 


06, 


.8 


0.59 


-1.0 


2, 


.6 


0.82 


N 


CH3OH-34 


20 


39 


01, 


.18 


42 


22 


06, 


.8 


0.88 


-3.7 


3, 


.3 


1.41 


N 


CH3OH-35 


20 


39 


01, 


.22 


42 


22 


41, 


.2 


0.58 


-3.0 


2, 


.6 


0.75 


Y-10 


CH3OH-36 


20 


39 


01, 


.47 


42 


23 


06, 


.3 


0.064 


-2.3 


1, 


.3 


0.073 


N, 


CH3OH-37 


20 


39 


01, 


.48 


42 


22 


07. 


.6 


0.78 


-4.3 


2, 


.0 


0.96 


N 


CH3OH-38 


20 


39 


01 


.50 


42 


22 


45. 


.2 


0.66 


-3.0 


2 


.0 


0.85 


Y-11 


CH3OH-39 


20 


39 


01, 


.51 


42 


22 


40. 


.8 


4.68 


-5.0 


3, 


.3 


6.07 


Y-12 


CH3OH-4O 


20 


39 


01, 


.54 


42 


22 


02, 


.8 


0.76 


-3.7 


2, 


.0 


1.00 


N 


CH3OH-4I 


20 


39 


01, 


.67 


42 


22 


43, 


.4 


2.50 


-6.3 


2, 


.6 


3.05 


Y-13 


CH3OH-42 


20 


39 


01, 


.70 


42 


22 


45. 


.2 


0.30 


-5.0 


1, 


.3 


0.38 


N 


CH3OH-43 


20 


39 


01, 


.76 


42 


22 


45. 


.3 


0.21 


-3.7 


1, 


.3 


0.20 


N, T 


CH3OH-44 


20 


39 


02, 


.02 


42 


22 


41. 


.2 


2.23 


-4.3 


1, 


.3 


2.62 


Y-14, 


CH3OH-45 


20 


39 


02, 


.02 


42 


22 


40. 


.8 


1.94 


-5.6 


1, 


.3 


1.82 


N, 
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Table 4 — Continued 



Maser ID 


RA (J2000) 


Dec (J2000) 


I" 

(Jyb-i) 


(km s-i) 


(km s-l) 


(Jyb-1 km s-l) 


Notes'' 


CH3OH-46 


20 39 02.08 


42 22 45.2 


6.29 


-3.7 


2.6 


6.45 


Y-15 


CH3OH-47 


20 39 02.21 


42 22 43.4 


0.32 


-4.3 


2.0 


0.30 


N, 


CH3OH-48 


20 39 02.22 


42 22 42.6 


1.82 


-3.7 


2.0 


1.80 


N, 


CH3OH-49 


20 39 02.24 


42 22 44.0 


8.77 


-5.0 


2.6 


12.5 


Y-17, 



^The typical rms in one channel is Smjybeam ^ for those channels with no strong (> 1 Jy) masers. The values 
reported in this table have been corrected for primary beam attenuation. AVso- is the full lincwidth at 3a level. 

''We mark with 'Y' or 'N' whether the maser was previously reported by Kurtz et al. (2004). The number after 
'Y' is the maser number as reported by Kurtz et al. (2004). We denote with 'O' maser features that axe blended in 
frequency with other maser emission. Tentative detections are marked with 'T'. 
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Table 5. Spectral Index between X and K Bands 





MMl-NW 


0.8±0.2 


MMl-SE 


0.8±0.2 


DR21(OH)-Rl 


-0.3±0.2 


DR21{OH)-R2 


0.1±0.2 


DR21(OH)-R4 


0.4±0.3 


DR21(OH)-R5 


-0.6±0.2 


DR21(OH)-R6 


0.4±0.2 



Notes: Spectral in- 

dex {Si, cc u") between 3.6 
and 1.3 cm using the values 
reported in Table 3. 
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01.25 0r.15 01.05 00.95 00.85 



20 39 01.4 01.3 01.2 01.1 01.0 00.9 00.8 00.7 00.6 00.5 

RIGHT ASCENSION (J2000) 

3.6 cm: Peak Intensity = 1.21 mJy/b 

Levs = 19 X (-3, 3, 5, 7, 10, 20, 30, 40, 50, 60, 63) nJy/b 
1.3 cm: Peak Intensity = 2.87 mJy/b 

Levs = 36 X (-3, 3, 4, 5, 10, 20, 30, 40, 50, 60, 70, 75, 78) [uiylb 
7 mm: Peak Intensity = 9.60 mJy/b 

Levs = 0.4 X (-3, 3, 5, 7, 9, 12, 16, 20, 22, 24) mJy/b 



Fig. 1. — Radio continuum emission in DR21(0H) at 3.6 cm (X-band), 1.3 cm (K-band) and 
7 mm (Q-band). The synthesized beams are shown in the lower left corner of the maps (see 
Table 2). The location of the 6.7 GHz CH3OH masers from Harvey-Smith et al. (2008) are 
shown with circles in the 1.3 cm continuum map; H2O masers are shown with diamonds in 
the 3.6 cm continuum map (Palmer & Goss, in prep.). 
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1401= 




5 
LSR Velocity (km/s) 



Fig. 2. — CH3OH 44 GHz spectrum of the DR21(0H) region. The spectrum was obtained 
from a 35" X 11" area encompassing the CH3OH maser arcs shown in Figure 4. 
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Fig. 3. — Integrated intensity (zero velocity moment) of the CH3OH 44 GHz masers in the 
DR21(0H) region. The numbers mark the position of the masers as given in Table 4. 
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Colors: Velocity Field of Methanol Masers 
Color Scale = [-7.0, 1 .0] km/s 
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02.0 



01.5 



01.0 00.5 00.0 
RIGHT ASCENSION (J2000) 

Contours: 3.6cm Continuum 

Contour Levels = 0.019 x (4, 10, 30, 60) mJy/b 



38 59.5 



59.0 



Fig. 4. — Velocity field (first velocity moment) of the CH3OH 44 GHz masers in the 
DR21(0H) region. The 3.6 cm continuum from Figure 1 is shown in contours; sources Rl 
and R6 are explicitly marked to facilitate comparison with Figure 1. The location of the 
NH3 molecular cores from Mangum et al. (1992) are shown with circles; the diameter of the 
circles equals the major axis of the cores as reported by Mangum et al. (1992). 
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10* 4 10* 

Freq (GHz) 



Fig. 5. — Radio spectral energy distribution of the MMl continuum detected in this work. 
The flux densities at X and K bands were obtained after convolving the maps to the synthe- 
sized beam of the Q band observations. A weighted power law fit {Si, oc u") to the data is 
shown. 
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Fig. 6. — Spitzer IRAC image of the DR21(0H) region (3.6 /xm blue; 4.5 /iin 
green; 8.0 /im red). The images were retrieved from the Spitzer archive 
(http:/ /irsa.ipac.caltech.edu/apphcations/Spitzer/Spitzer/ ). The crosses mark the positions 
of CH3OH 44 GHz masers from this work (Figure 3) and the two circles show the positions 
of the radio sources MMl-NW and SE (Figure 1). Note the 8.0//m dark area oriented north- 
south that corresponds to the large scale molecular filament in the DR21 region (e.g., Kumar 
et al. 2007). 



